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Abstract We give an overview of the experimental measurements and the theoretical un- 
derstanding of the branching fractions and CP-violating asymmetries of charmless B-meson 
decays. Most experimetal results are from the BABAR and Belle experiments during the past 
decade. The global features of these experimental results are typically well described by the 
QCD-motivated theories such as QCD factorization, pQCD and soft-collinear effective theory. 
The agreement between theory and experiment is generally satisfactory, though there remain 
some unsolved puzzles that pose a great challenge to both theorists and experimentalists. 
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1 INTRODUCTION 

Evidence for the B meson was first seen in 1981 (1). For the next two decades, the 
ARGUS experiment operating at DESY in Hamburg, Germany and the CLEO 
experiment operating at CESR at Cornell University studied the properties of 
these mesons and made many important discoveries including the discovery of 
B°B° mixing. The first examples of charmless hadronic B decays were seen by 
CLEO (2). Starting in 2000, there have been many measurements of these decays 
with increasing precision by the BABAR experiment operating at PEP-II at 
SLAC in California and the Belle experiment operating at KEK in Japan. There 
are now nearly 100 of these decays that have been observed with a statistical 
significance of at least four standard deviations. While the early measurements 
by CLEO were groundbreaking, the errors are sufficiently large that they have 
little weight in the present world averages. The CDF experiment at FNAL also 
has measurements for these decays. While they are nearly as precise as BABAR 
and Belle, there are only measurements for four decay channels. Thus in tables 
throughout this review, we concentrate on the copious measurements from Belle 
and BABAR. 

The theoretical study of weak nonleptonic decays of a heavy meson is diffi- 
cult and involved due to the interplay of short- and long-distance QCD effects. 
The effective weak Hamiltonian at the quark level H = c i(^)Oi(/J-) is theoret- 
ically well under control, where Oi are four-quark operators and Ci(/J,) are the 
Wilson coefficients which incorporate strong-interaction effects above the scale 
fj,. However, it is a difficult task to evaluate the hadronic matrix elements of the 
local operator Oi reliably due to the nonperturbative QCD effects involved. A 
simple and widely employed approach is based on the valence quark assumption 
and the vacuum-insertion (or factorization) approximation in which the hadronic 
matrix elements of two quark bilinear operators are saturated by the vacuum 
intermediate states. 

Due to the experimental and theoretical efforts in the past decade, qualitative 
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understanding of nonleptonic charmless B decays has become possible. Since 
the B meson is heavy, it is possible to describe the dynamics of hadronic decays 
by theories motivated by QCD rather than by phenomenological models. A 
central aspect of those theories is the factorization theorem which allows us to 
disentangle short-distance QCD dynamics from nonperturbative hadronic effects. 
In the heavy quark limit, matrix elements can be expressed in terms of certain 
nonperturbative input quantities such as light cone distribution amplitudes and 
transition form factors. Power corrections beyond the heavy quark limit generally 
give the major theoretical uncertainties. 

In this article we give an overview of the experimental measurements and the 
theoretical understanding of the branching fractions and CP-violating asymme- 
tries of charmless P-meson decays. We begin with the theoretical and experi- 
mental tools necessary for the study of hadronic B decays. This is followed with 
a discussion of 2-body, 3-body, quasi-2-body, and baryonic B decays, and the sta- 
tus of time-dependent CP-violation measurements and predictions. We discuss 
the puzzles that remain unsolved. 

2 FACTORIZATION AND THEORETICAL TOOLS 

In this section we will introduce various approaches that have been employed 
for studying the dynamics of hadronic B — > M\M 2 decays. In the effective 
Hamiltonian approach, the decay amplitude is given by 

A(B - MxM 2 ) = ^ E AiCi^XMiMalOilBX/i), (1) 

where Aj are Cabibbo-Kobayashi-Maskawa (CKM) (3) matrix elements, Oi are 
four-quark operators and Cj(/x) are the Wilson coefficients which incorporate 
strong-interaction effects above the scale /x. A major theoretical issue is how 
to evaluate the matrix elements of the four-quark operators (M\M2\Oi\B). 

2.1 Naive and Generalized Factorization 

A widely used approximation is the so-called "naive factorization" or "vacuum- 
insertion approximation" under which the matrix element [M\M 2 \0\B) is ap- 
proximated by {M 1 \J 1 ^\0){M 2 \J^\B) or (M 2 \Ji tl \0)(M 1 \J^\B) with being a 
bilinear current; that is, the matrix element of a four-quark operator is expressed 
as a product of a decay constant and a form factor. Naive factorization is sim- 
ple but fails to describe color-suppressed modes. For example, the predicted 
ratio of T(D° -► K°tt )/T(D° -► K~tt + ) « 3 x 10" 4 is too small compared 
with the experimental value of 0.55 (4). This is ascribed to the fact that color- 
suppressed decays receive sizable nonfactorizable contributions that have been 
neglected in naive factorization. Another issue is that the decay amplitude under 
naive factorization is not truly physical because the renormalization scale and 
scheme dependence of Cj(/z) are not compensated by that of the matrix element 
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(MiM2\Oi\B)(fj,). In the improved "generalized factorization" approach (5,6), 
nonfactorizable effects are absorbed into the parameter N£ s , the effective num- 
ber of colors. This parameter can be empirically determined from experiment. 
Since these early calculations have been replaced with improved ones discussed 
below, we will not discuss these early calculations. 

2.2 Theories of Hadronic B Decays 

With the advent of heavy quark effective theory, nonleptonic B decays can be 
analyzed systematically within the QCD framework. There are three popular ap- 
proaches available in this regard: QCD factorization (QCDF) (7), perturbative 
QCD (pQCD) (8) and soft-collinear effective theory (SCET) (9). A detailed dis- 
cussion of these theories goes beyond the scope of this review, and the interested 
reader is referred to the original literature. Basically, theories of hadronic B de- 
cays are based on the "factorization theorem" under which the short-distance con- 
tributions to the decay amplitudes can be separated from the process-independent 
long-distance parts. 1 In the QCDF approach, nonfactorizable contributions to 
the hadronic matrix elements can be absorbed into the effective parameters 

A(B - M X M 2 ) = ^J2 X MMiM 2 )(M 1 M 2 \O l \B) iact , (2) 

where cij are basically the Wilson coefficients in conjunction with short-distance 
nonfactorizable corrections such as vertex, penguin corrections and hard spectator 
interactions, and (MiM2\Oi\B)f a , ct is the matrix element evaluated under the 
factorization (or vacuum insertion) approximation. For penguin operators Oq$, 
the current operators are replaced by the scalar or pseudoscalar densities. 
Since power corrections of order Aqcd/^6 are suppressed in the heavy quark 
limit, nonfactorizable corrections to nonleptonic decays are calculable. In the 
limits of m& — ► oo and a s — > 0, naive factorization is recovered in both QCDF 
and pQCD approaches. 

Power corrections are often plagued by the end-point divergence that in turn 
breaks the factorization theorem. For example, the endpoint divergence occurs 
at the twist-3 level for the hard spectator scattering amplitude and at the twist-2 
level for the annihilation amplitude. As a consequence, the estimate of power cor- 
rections is generally model dependent and can only be studied in a phenomeno- 
logical way. There is also an endpoint singularity in SCET though this issue 
can be resolved after introducing the zero-bin subtraction procedure (10). In 
the pQCD approach, the endpoint singularity is cured by including the parton's 
transverse momentum. 

2.3 Diagrammatic Approach 

Because a reliable evaluation of hadronic matrix elements is very difficult in 
general, an alternative approach is based on the diagrammatic approach. It has 

1 However, the charming penguin term advocated in SCET violates this factorization theorem. 
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been established that a least model-dependent analysis of heavy meson decays 
can be carried out in the so-called quark-diagram approach (11, 12). In the 
diagrammatic approach, all two-body nonleptonic weak decays of heavy mesons 
can be expressed in terms of six distinct quark diagrams (Fig. I): 2 T, the color- 
allowed external W-emission tree diagram; C, the color-suppressed internal W- 
emission diagram; E, the W-exchange diagram; A, the W-annihilation diagram; 
P, the penguin diagram; and V, the vertical TV-loop diagram. It should be 
stressed that these quark diagrams are classified according to the topologies of 
weak interactions with all strong interaction effects included and hence they are 
not Feynman graphs. All quark graphs used in this approach are topological with 
all the strong interactions included, i.e. gluon lines are included in all possible 
ways. The diagrammatic approach was applied to hadronic B decays first in (14). 
Various topological amplitudes have been extracted from the data in (13, 15-17) 
after making some reasonable approximations, e.g., SU(3) symmetry. 




T C E 




A P V 



Figure 1: Various topological diagrams for B — ► M1M2 decays. 



3 EXPERIMENTAL TOOLS 

The experimental measurements involve separation of small samples (10 — 2000 
signal events) from total samples of BB and light quark-pairs (qq) of several 
billion events. The background is typically dominated by the copious qq produc- 
tion, where the background after preliminary sample selection is often 1000 times 

2 Historically, the quark-graph amplitudes T, C, E, A, P named in (13) were originally de- 
noted by A, B, C, D, E, respectively, in (11, 12). For the analysis of charmless B decays, one 
adds the variants of the penguin diagram such as the electroweak penguin and the penguin 
annihilation. 
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larger than the signal. The separation of the signal from the large background 
relies on a number of variables that are common to all analyses. For signal, the 
mass of a candidate, reconstructed from the charged and neutral tracks in the 
event, is equal to the B mass, approximately 5.28 GeV/c 2 . The center-of-mass 
energy of the candidate is equal to one-half of the T(45) rest energy since all e + e~ 
B experiments take advantage of the increased production rate of BB events at 
this resonance. Finally there are several variables that take advantage of the dif- 
ference in the shape of the signal and background events. In the center-of-mass 
system, signal is typically spherical since the B mesons are produced nearly at 
rest. The qq background events are are characterized by back-to-back "jets" of 
particles with small transverse momentum with respect to the direction of the 
leading quarks. For B decays involving resonances in the final state (p, K*, 77, rj', 
u>, 4>, etc), the mass of the daughter particles from the decay of these resonances 
is also used. This is useful because the background often is dominated by events 
with combinatorial background, where the resonance candidate is not real but 
rather is composed of combinations of particles from the qq event that happen to 
have an invariant mass near that of the resonance. For spin-1 (vector) particles, 
the so-called helicity angle is also often useful since the decay of the daughters in 
the vector particle's rest frame is not typically uniform. 

Most analyses combine some of the above quantities into a maximum likelihood 
(ML) fit. Such a fit characterizes the signal and background with probability den- 
sity functions (PDFs) that describe the distribution expected for each variable. 
It is the difference between the shapes of the signal and background PDFs that 
allow extraction of the signal from the very large backgrounds. The shape of 
these PDFs is typically determined from Monte Carlo (MC) simulation for signal 
and BB background events, with checks provided by "control" samples of other 
more copious decays. The data itself is used to determine the PDFs for the qq 
background. The free variables in the ML fit typically include the yields of the 
signal and backgrounds, the value of CP violation parameters where relevant, 
the longitudinal polarization fraction for decays with two particles with non-zero 
spin, and often some of the parameters that determine the qq background shapes. 
The structure of the likelihood often assumes that the input observables are un- 
corrected. This assumption is tested with data and MC and correlations are 
typically below 10%. The residual correlations may cause small signal biases 
(~10%) which are evaluated with MC and appropriate corrections are made. 

4 TWO-BODY DECAYS 
4.1 Branching Fractions 

The general expressions of topological amplitudes are 
A(B°^7r + 7r-) = T + P+^P£ W + E + V, 
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A(B°-7r 7r°) = ~(C - P + P EW + - E - V), (3) 

- tt+tt ) = -L(T + C + P EW + P£ W ), 

for tree-dominated B — > tttt decays, 
A(B° -> K + K~) = E + P A , 
A(B ^K K°) = P-lp^ + P A , 

A(B+^K+K°) = A + P-±P£ W , (4) 
for B — > Kif decays, and 

A(B°-ffV) = -^(P'-C'-P^w-^w + ^a), (5) 

A(P + ^K°7T + ) = P'-lp^ w +A' + PA, 

A(B + ^K+n°) = -L(P> + T' + C' + P EW + lpg w + A' + P A ), 

for B — ► i^7r decays, where Pew an d Pew are color- allowed and color-suppressed 
electroweak penguin amplitudes, respectively, and Pa is the penguin-induced 
weak annihilation amplitude. We use unprimed and primed symbols to denote 
AS = and \AS\ = 1 transitions. 

Experimental results for branching fractions from BABAR and Belle and the- 
oretical predictions are summarized in Fig. 2. Here and in the following tables, 
the theoretical values and errors are from weighted averages of the various pre- 
dictions with the errors divided by y/3 since the quoted theory errors indicate a 
range as parameters are varied, not la errors. 

As mentioned in Sec. 2.2, endpoint divergences will occur in the QCDF ap- 
proach in the penguin annihilation and hard spectator scattering amplitudes, 
which are often parametrized as (7) 

X A = In (1 + PA e^), X H = In (1 + p H e^), (6) 

where the parameters pa,h and 4>a,h are real and A^ is a typical hadronic scale of 
order 500 MeV. Since these parameters are unknown within QCD factorization, 
the central values of QCDF predictions correspond to pa,h = and the power 
corrections due to annihilation effects give the largest theoretical uncertainties. 
Hence, for penguin-dominated modes, the central values predicted by QCDF 
are usually systematically below the measurements. Penguin annihilation and 
hard spectator scattering amplitudes are calculable in the pQCD approach as 
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the endpoint singularity is overcome by the parton's transverse momentum. The 
predicted central values by pQCD for penguin-dominated decays are normally 
higher than those of QCDF. 

The predicted 7r + 7r~ branching fraction of ~ 7 x is too large, whereas 

of order 0.3 x 10 6 is too small compared with experiment (the ratio of 
the branching fractions is predicted much more precisely than one would infer 
from Fig. 2). If the penguin amplitudes P, Pew, Pa an d the annihilation E are 
neglected, a fit to the ratio of B° — > 7r°7r° and B — > 7r + 7r~ rates will yield \C/T\ ~ 
0.50 — 0.60. In the short-distance (SD) factorization approach, |C/T|sd ~ 0.2 — 
0.3. It is thus a challenge to theorists to understand B° —* 7r°7r° and 7r + 7r~ 
decays. 




Branching Ratio x 10 6 

Figure 2: Branching fraction measurements of 2-body decays from BABAR (blue) 
(18-21) and Belle (red) (22-24) and theoretical predictions (green) (25-28). 

As for B -» KK decays, B+ -> K+K° and B° -> K°K° are dominated by the 
CKM-suppressed b —> dg penguin, while B° — > proceeds only through 

weak annihilation. Hence, the branching fraction is expected to be of order 10~ 6 
for the former and 10 -8 for the latter, in agreement with the data. 

The B — ► Kit decays are dominated by penguin contributions because of 
^ l^cs^SI ~ l^ts^tbl and the large top quark mass. For the ratios 
defined by 

= 2T(B+ -+ K+^) = T(B°^K+7r-) 

Uc - T (B+ if°vr+) ' n - 2T{B* -> ^vrO) ' lU 

we have i? c = if the other quark-diagram amplitudes are negligible compared 
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with P' . The current experimental measurements give R c = 1.12 ± 0.07 and 
R n = 0.99 ± 0.07. There are two approximate sum rules for B — ► Kit rates (29) 
and rate asymmetries (30) 

t{k + tt-) + t(k tt + ) « 2[r(irV) + r(]fV)], 

Ar(^+vr-) + Ar(KV+) « 2[Ar(K + 7r°) + Ar( J PC°7r )], (8) 

based on isospin symmetry, where AT(Ktt) = F(B — > K7f) — r(S — ► K7r). The 
rate sum rule is fairly well satisfied by experiment. 



4.2 Direct CP Asymmetries 

It is well known that it requires at least two amplitudes with nontrivial rela- 
tive strong and weak phases to produce direct CP violation. The direct CP 
asymmetry observed in B° — ► K + ir~ decays, 

T(B° ^ K'tt+)+T(B° ^ K+TT-) 

requires that the relative strong phase S between t' = T' + -P^w ano - p' = 
P> ~ 1 P ew + P A t see Eq- ( 5 )] be of order 15°. In QCDF, the phase S induced 
perturbatively from vertex corrections and penguin diagrams is too small and 
has a wrong sign (27). Therefore, it is necessary to include the contribution from 
the long-distance (LD) strong phase. The nonperturbative strong phase induced 
from the weak decay B — ► D^D^ followed by the final-state rescattering of 
£)(*)_£)(*) — > i^7r can reproduce the experimental observation. In the pQCD ap- 
proach (8), a large SD strong phase arises from the penguin-annihilation diagram 
P' A , though this is in contrast to the conventional wisdom that strong phases are 
basically nonperturbative in nature. 

In Fig. 3 we show Acp for the most precisely measured decays, typically those 
with an average uncertainty <0.10 or a deviation from zero of at least three 
standard deviations. Of these only B° — > K + ir~ and B + — ► p°K + have world 
averages for Aqp that are more than 4 standard deviations from 0. Theory pre- 
dictions are not included since the lack of knowledge of strong phases makes 
them very uncertain. All of these decays are measured from charge asymmetries 
as suggested by Eq. 9. We defer until Sec. 8 a discussion of the CP asym- 
metries of CP-eigenstate modes such as B° — > 7r + 7r~ that are measured with 
time-dependent techniques. 

If the color-suppressed C, color-allowed electroweak penguin Pg W and annihi- 
lation A' are negligible compared with the tree amplitude T', it is clear from Eq. 
(5) that the decay amplitudes of K + tt° and K + tt~ will be the same apart from a 
trivial factor of l/y/2. Hence, one will expect that Acp(K + 7t°) ~ Aq-p(K + tt~), 
while they differ by 5.3 a experimentally, AAkw = Acp(K + tt°) — Acp(K + it~) = 
+0.148±0.028. Since A(B + -► i^T+vr ) oc t' + c'+p' and A(B° -► K + tt~) oc t'+p' 
with d = C' + P^ff, the puzzle is resolved provided that d /t' is of order 0.5 ~ 0.6 
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TT+TT 

K + ir + Tr~ 



K+TT- 
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+ 



p~K+ 



p°K- 



0.7 



CP Asymmetry 



Figure 3: A sample of the most precise direct CP measurements from BABAR 
(blue) and Belle (red). All of these measurements are obtained from time- 
independent charge asymmetries. 
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with a negative relative phase. There are two possibilities for a large d: either 
a large color suppressed C or a large electroweak penguin Pg W . We note that a 
global fit to 7T7T, Kit and KK data in the diagrammatic approach gives a relative 
strong phase of (— 56± 10)° between C' and T' (16). Various scenarios for accom- 
modating large C (28) or (31) have been proposed. It is evident from the 
discussion above that a large color-suppressed amplitude with a sizable relative 
phase to the tree amplitude can solve both the Kir puzzle in CP violation and 
the enhancement of B° — > 7r°7r° in rate. 

Based on SU(3) flavor symmetry, direct CP asymmetries in Kn and tttt systems 
are related as (32): 

Ar^+o = -Ar^o, Ar^V) = -Ar^V). (10) 

The first relation leads to A C p(ir + Tr-) = [B(K + ir~) /B(Tr + ir-)]A C p(K+ir-) « 
0.37 , which is in good agreement with the current world average of 0.38 ± 0.06. 



5 THREE-BODY DECAYS 



Three-body decays of heavy mesons are more complicated than the two-body 
case as they receive both resonant and nonresonant contributions. They are 
generally dominated by intermediate vector and scalar resonances. The analysis 
of these decays using the Dalitz plot technique enables one to study the properties 
of various resonances. Moreover, the Dalitz plot analysis of 3-body B decays 
provides a nice methodology for extracting information on the unitarity triangle 
in the standard model. The experimental results are summarized in Fig. 4. The 
agreement between Belle, BABAR, and theoretical predictions is good in all cases. 



5.1 Resonant Contributions 

Under the factorization approximation, vector and scalar meson resonances con- 
tribute to the two-body matrix elements {P\P 2 \ V^|0) and (P\P2\S\Q), respectively, 
with S being a scalar density. They can also contribute to the three-body matrix 
element {PiP^V^ — A^B). Resonant effects are generally described in terms of 
the usual Breit-Wigner formalism. For example, 

(K + ( Pl )K-(p 2 )\q^ q \0) R = ^(K+K-lVi) — ( Vi \q^q\0), 

i m Vi ~ s ~ lm Vi L Vi 

(K + ( Pl )K-(p 2 )\ss\0) R = J2(K + K-\S t ) — — (SiMO), (11) 

i ™ s .- s-im Si Tsi 

where s = ( Pl +p 2 ) 2 , Vi = </>,p,u, ■ ■ ■ and S t = / (980), / (1370), / (1500), • ■ ■. 
Once the strong couplings for Vi(Si) — > K + K~ and the decay constants of Vi 
and Si are known, we are ready to compute various resonant contributions to the 
3-body decays of interest. Using the narrow width approximation for resonance 
R 

F(B -^RP -> P!P 2 P) = T(B -► RP)B{R -► PiP 2 ), (12) 
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Figure 4: Branching fraction measurements of 3-body decays from BABAR (blue) 
(33-44) and Belle (red) (45-48) and theoretical predictions (green) (49). The 
predictions include both resonant and nonresonant contributions. 

we can calculate the rates for the quasi- two-body decays B — > PV and B — > SP. 
For the details of theoretical calculations of charmless 3-body decays of B mesons 
based on the factorization approach, see (49); for a theoretical overview, see (50). 

5.2 Nonresonant Contributions 

One of the salient features of 3-body decays is the large nonresonant contribution 
to penguin-dominated modes. It is known that the nonresonant signal in charm 
decays is small, less than 10% (51). On the contrary, the nonresonant fraction is 
about ~ 90% in B — > KKK decays, ~ (17—40)% in B — > Ktttt decays (smaller in 
the Ktttt decay), and ~ 14% in the B — > tttttt decay. Hence, the nonresonant 3- 
body decays play an essential role in penguin-dominated B decays. Nonresonant 
amplitudes in charm decays are usually assumed to be uniform in phase space. 
However, this is no longer true in B decays due to the large energy release in 
weak B decays. This makes the Dalitz plot analysis of nonresonant contributions 
more difficult. While both BABAR (37, 41) and Belle (46) have adopted the 
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parametrization 

^nr = (c 12 e^ 12 e- as ^ + ci 3 e^ 13 e" QS ?3 + c 2 3e^ 23 e~ as ^)(l + 6 N Re i(/3+<5NR) ) (13) 

to describe the non-resonant B — > KKK amplitudes, they differ in the analysis 
of the nonresonant component in B — ► i^7T7r decays. 

A detailed theoretical analysis in (49) indicates two distinct sources of nonres- 
onant contributions: a small contribution from the tree transition and a large one 
from the matrix elements of scalar densities, e.g., (KK\ss\0), induced from the 
penguin transition. This explains the dominance of the nonresonant background 
in B — > KKK decays, the sizable nonresonant fraction in K~ir + ir~ and K ir + ir~ 
modes and the smallness of nonresonant rates in B — > tttttt decays. 

Mixing-induced CP asymmetries of 3-body decays will be discussed in Sec. 7 
below. 

6 QUASI-2-BODY DECAYS 
6.1 B -> PP 

Among the 2-body B decays, B — > r/.fC has the largest branching fraction, of 
order 70 x 1CT 6 , while B(B — > 77-ff) is only (1 - 3) x 10~ 6 . This can be roughly 
understood as follows. Let's express the rf and rj wave functions in the quark- 
flavor basis r\ q = (uu + dd) / \/2 and rj s = ss 

r] = cos (j)rj q — sin cpr] s , rf = sin 4>rj q + cos ^ry s , (14) 

where the mixing angle is extracted from data to be ^ = 39.3° (52). The in- 
terference between the B — > amplitude induced by the 6 — > sgg penguin 
and the i? — > ?? s i^ amplitude induced by 6 — > sss is constructive for i? — ► r/'i^ 
and destructive for B — > 77^. This explains the large rate of the former and the 
suppression of the latter. However, the predicted rates are still somewhat smaller 
than the measurements (see Fig. 5). 

Many possible solutions to the puzzle for the abnormally large rj'K rate have 
been proposed: (i) a significant flavor-singlet contribution (15,16,64), (ii) a large 
B —> 7]' form factor (68), (iii) a contribution from the charm content of the rf ', (iv) 
an enhanced hadronic matrix element (0|s75s|r/) due to the axial U(l) anomaly 
(69), (v) a large chiral scale rriQ associated with the r] q (70), (vi) a long-distance 
charming penguin in SCET (26), and (vii) a large contribution from the two-gluon 
fusion mechanism (71). 

Because the rj' is dominated by the flavor-singlet component 771, it is plausible 
that a sizable flavor-singlet amplitude S' will account for the large r/'K rates. 
When this contribution is included in the diagrammatic approach, a global fit 
to the B — > PP data leads to a flavor-singlet amplitude S about 3 to 4 times 
the magnitude of -P^w (15,16). This flavor-singlet contribution is also supported 
by the consideration of the B — > r/'V° decays. Their topological quark decay 
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Figure 5: Branching fraction measurements of decays with rj and rj' mesons 
from BABAR (blue) (53-58) and Belle (red) (59-63) and theoretical predictions 
(green) (26,27,64-67). 



amplitudes are 

A(B + 777T+) = — ]=(t + c + 2p + s), A(B° -> rjTT ) = — l —(2p - s), 
V3 v 6 

A(B + -► 7?'^+) = -i=(t + c + 2p + 4s), ^(B -► t/tt ) = -L(p + 2a), (15) 

where i = T + P|w c = C + P E w, p = P - §-Few + Pa, s = S - |ffew, and 
we have assumed = 35.3° since this value is algebraically simple. Note that 
there is no tree contribution to the r/V modes. If s = 0, one will naively 
expect that B{B° -» t/tt ) = \B(B° r?7r°). QCDF predicts 5(5° r/vr°) = 
(0.28±^;||) x 10" 6 and 8(5° -► t/tt ) = (0.17±g;ff) x 10~ 6 (27). The pQCD 
approach (65) has very similar results. From the experimental measurements 
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B(B° -► r/vr°) < 1.5 x 10" 6 and B(B° -► t/tt ) = (0.9±0.4) x 10~ 6 by BABAR (56) 
and 6(5° -► t/tt ) = (2.8 ± 1.0) x 10~ 6 by Belle (60), it appears that the current 
predictions of QCDF and pQCD may be too small for B° — > 7/71" . In the presence 
of S, B° -► t/tt is enhanced: 6(5° -► 77V ) « ^(5° -► r/vr ) ~ 1.0 x 10~ 6 (16). 

Since the two penguin processes b — > sss and 6 — > sgg contribute destructively 
to £? ^ ry-RT, the penguin amplitude is comparable in magnitude to the tree ampli- 
tude induced from b — > usu, contrary to the decay S — > 7/ if which is dominated 
by large penguin amplitudes. Consequently, a sizable direct CP asymmetry is 
expected in B + — > r?K + but not in r]'K + (72). Indeed, the average of BABAR 
and Belle measurements yields a 3<r effect, Acp{rjK + ) = —0.27 ± 0.09 . 

6.2 B -> PV 

The experimental results are summarized in Fig. 6. The decays B — > i-V have 
been studied in QCD-motivated approaches such as QCDF, pQCD and SCET; 
see (66) for comparison of the theory predictions in various approaches with 
experiment. For tree dominated decays, the predictions of QCDF and pQCD are 
in agreement with experiment for B° — > /9 ± 7r =F and B — ► wzr, while the SCET 
calculations are smaller due to the smallness of both B — > 7r and i? — > p form 
factors in SCET. On the contrary, SCET predicts larger rates than QCDF and 
pQCD for the color-suppressed decays such as B° — > p°7r°, rf^fp, rf ^uj because 
the hard form factor £j 17 is comparable with the soft part £ (F PV = Qj V + 
C Py ) and is enhanced by a large Wilson coefficient. The predicted B(B° — > /9°7r°), 
of order 0.4 x 10~ 6 by QCDF (27) and of order 0.1 x 10~ 6 by pQCD (90), are 
too small compared with the experimental average of (2.0 ± 0.5) x 10 -6 while 
the SCET calculations which rely on some input from experiment are consistent 
with experiment (66). 

For penguin-dominated B — > PV decays, the predictions by QCDF are sys- 
tematically below the measurements (27), while the data can be accommodated 
by fitting them to SCET. 

As for A C p(B° -»• K + it~), QCDF predicts a wrong sign for A C p(B + -► p°K+). 
The pQCD approach predicts too large direct CP violation in many PV modes 
(66). 

6.3 B VV 

The underlying mechanism for B — > VV is more complicated than PV and 
PP modes as it involves three polarization vectors. The decay amplitude of 
B — > yy can be decomposed into three components, one for each helicity of the 
final state: Aq,A + ,A-. The transverse amplitudes defined in the transversity 
basis are related to the helicity ones via 



A 



A+ + A 
V2 



A+ - A_ 



V2 



(16) 
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Figure 6: Branching fraction measurements of PV decays from BABAR (blue) 
(35,37,43,44,56,57,73-82) and Belle (red) (46-48,83-89) and theoretical predic- 
tions (green) (27,66,90-93). 



It is naively expected that the helicity amplitudes Ah (h = 0, +, — ) for B — > VV 
respect the hierarchy pattern Aq : A + : A- = 1 : (Aqcd/^6) : (Aqcd/^) 2 - 
Hence, charmless B — > VV decays are expected to be dominated by the longitu- 
dinal polarization states and satisfy the scaling law, 



i-f L =o { m, i±=i +0 (^), d7) 

\m 2 B J f\\ \m B J 
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Figure 7: Longitudinal polarization measurements from BABAR (blue) (99-109) 
and Belle (red) (84,110-113) and theoretical predictions (green) (114-117). 



with fi, fx and fn being the longitudinal, perpendicular, and parallel polariza- 
tion fractions, respectively, defined by 

Y a \A a \ 2 
fa = T = |Ao| 2 + |A||| 2 + |Al| 2 ' (18) 

with a = L, ||, -L. The experimental measurements are summarized in Fig. 7. 
The progress for these measurements is quite impressive, with the uncertainty 
for a handful of modes now < 0.05 . In sharp contrast to the pp case, the large 
fraction of transverse polarization observed in B — > K*p and B — > <pK* decays 
at B factories is a surprise and poses an interesting challenge for theoretical 
interpretations. Various mechanisms such as sizable penguin-induced annihilation 
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contributions (94), 3 final-state interactions (95,96), form-factor tuning (97) and 
new physics (98) have been proposed for solving the B — ► VV polarization puzzle. 

Two recent calculations (116,117) indicate that NLO nonfactorizable correc- 
tions from vertex and penguin corrections and hard spectator scattering will 
render the positive-helicity amplitude of some VV modes comparable to the lon- 
gitudinal one and hence will increase the transverse polarization. For example, 
/l is naively expected to be 1 — Amy /m 2 B ~ 0.90 in B — > tfiK* and B° — > K*°p° 
decays. However, NLO corrections decrease this expectation to /^((pK*) ~ 0.6 
and fL(K*°p°) ~ 0.5. Therefore, the polarization puzzle is alleviated to a large 
extent by the consideration of NLO effects. The theoretical predictions in Fig. 7 
reflect these more recent calculations. 

According to the recent calculations based on QCDF (116, 117), there is a 
hierarchy pattern for the polarization fractions in B — ► K*p decays: 

f L (K* + p°) > f L (K* + p-) > f L (K*°p+) > f L (K*°p°). (19) 

This pattern is compatible with measurements though only two, K*°p + and 
K*°p°, are well measured. Improved measurements of all of these decays are 
important in further testing the theoretical calculations. 

Even though Jl can be substantially reduced in the presence of nonfactoriz- 
able corrections, the polarization anomaly is not fully resolved unless the rate 
is also reproduced correctly. The experimental branching fraction measurements 
are summarized in Fig. 8. In most cases the agreement between theory and ex- 
periment is quite good. The recent measurements of the B — ► toK* decays are 
well below the predicted average though this is somewhat misleading since the 
average is mostly from pQCD (115); the QCDF prediction (116,117) is about a 
factor of two smaller and in reasonable agreement with the data. The B — > cpK* 
rate predicted by QCDF (pQCD) is too small (large) compared with the data. To 
improve the situation, QCD factorization and pQCD (121) rely on penguin an- 
nihilation amplitudes, while SCET invokes charming penguins (9), and the final- 
state interaction model considers final-state rescattering of intermediate charm 
states (95,96). 

In QCDF, the theoretical model predicts B(B° p°p°) ~ 0.9 x 10~ 6 (116,117) 
and B(B° — > 7r°7r°) ~ 0.3 x 10~ 6 (27), whereas experimentally the latter has a rate 
larger than the former. One plausible possibility is that final-state interactions 
are important for B — ► -ktt but not for B — > pp. The B — > tt-k amplitudes can 
be decomposed into the 1 = and 2 isospin states with isospin phases <5g and 
5% , respectively. When the isospin phase difference is sizable, the 7r°7r° mode will 
be enhanced by the final-state rescattering of 7r + 7r - to 7r°7r° (this amounts to 
enhancing the color-suppressed amplitude C, see also (123)). Since B — > p + p~ 
has a rate much larger than B — ► it + it~, it is natural to expect that B — > p°p° 
will receive large enhancement via isospin final-state interactions. The fact that 
the branching fraction of this mode is rather small and is consistent with the 

3 Historically, even before the experimental observation of the polarization puzzle, it was 
already pointed out in (118) that penguin annihilation effects could reduce fL(4>K") to 0.75. 
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Figure 8: Branching fraction measurements of VV decays from BABAR (blue) 
(100-106,108,109) and Belle (red) (84,110-113,120) and theoretical predictions 
(green) (114-117,122). 



theory prediction implies that the isospin phase difference of 5q and 5% must be 
negligible and so is the final-state interaction (124). 

6.4 B -> (5, A, T)M 

Much less is known about scalar, axial-vector and tensor mesons, though there 
are already dozens of measurements involving these modes as decay products of 
B mesons. In this section we summarize the experimental measurements (see 
Fig. 9) and the theoretical expectations. 

6.4.1 B — ► SP, SV It is known that the identification of light scalar 
mesons is difficult experimentally and the underlying structure of scalar mesons is 
not well established theoretically. It is hoped that through the study of B — > SP, 
old puzzles related to the internal structure and related parameters, e.g. the 
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Figure 9: Branching fraction measurements of decays involving + (bottom 
group), 1 + (middle group), and 2 + mesons (top group) from BABAR (blue) 
(36,43,44,53,101,107,108,126-131) and Belle (red) (46-48,132) and theoretical 
predictions (green) (117,125,134-139). 
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masses and widths, of light scalar mesons can receive new understanding. For ex- 
ample, it has been shown (125) that if ao(980) is a qq bound state, the predictions 
are B(B° -► a^(980)vr^) ~ 8.2 x 1(T 6 and B(B° -► a (980)K+) ~ 4.3 x lO" 6 . 
These exceed the current experimental 90% confidence level (C.L.) upper limits 
3.1 x 1CT 6 and 1.9 x 10 -6 (54), respectively, suggesting that the four-quark nature 
for the ao (980) is favored. 

One of the salient features of the scalar meson is that its vector decay constant 
fs defined by (S(p)|g27/i9i|0) = fsP/j, is either zero or small (of order — m u , 
m s — m^u). Therefore, when one of the pseudoscalar mesons in B — > PP decays 
is replaced by the corresponding scalar, the resulting decay pattern could be 
very different. For example, it is expected that T(B + — > ao~7r°) <C T(B + — > 
ag7r + ) and T(B° — > (Iq1t~) <C T(B° — > Oq 7r + ) as the factorizable contribution 
proportional to the decay constant of the scalar meson is suppressed relative to 
the one proportional to the pseudoscalar meson decay constant. This feature can 
be checked experimentally. 

The decay B — ► /o(980)i^T is the first charmless B decay into a scalar meson 
observed at B factories (133). It receives two different types of penguin con- 
tributions: one from b — ► suu and the other from b — > sss. Due to the large 
scalar decay constant /| Q of order 370 MeV defined by (/o|ss|0) = fnf ff that 
appears in the penguin amplitude, this decay is dominated by the b — ► sss pen- 
guin contribution with predictions of about 15 x 10~ 6 (125). The experimental 
measurements are for B(B — > / (980)if) x B(/ o (980) — > 7r + 7r~). This product is 
of order (5 — 10) x 10 -6 for the B + and B° decays. The theoretical predictions 
are consistent with experiment provided that #(/o(980) — > 7r + 7r~) ~ 0.50. 

6.4.2 S — ► AP There are two distinct types of parity-even axial- vector 
mesons, namely, 3 Pi and l P\. The 3 -Pi nonet consists of ai(1260), /i(1285), 
/i(1420) and K 1A , while the l P x nonet has &i(1235), ^i(1170), /ii(1380) and 
K\B- The physical mass eigenstates if i (1270) and -fTi(1400) are mixtures of K\a 
and K\b states owing to the mass difference of the strange and non-strange light 
quarks. 

A prominent feature of the 1 P\ axial vector meson is that its axial- vector decay 
constant is small, vanishing in the SU(3) limit. This feature was confirmed by the 
BABAR observation (127) that T(B° -► 6+vr-) <C T(B° -► 67/ tt + ). By contrast, 
it is expected that T(B° -> a^ - ) > r(B° -> af 7r+) since / ai > /„-. 

The predicted branching fractions for the b\K and b\Ti modes are in good agree- 
ment with the BABAR measurements. The comparison for of ir^ 1 and a^K + is 
sometimes good and sometimes not so good (see Fig. 9); improved measurements 
are needed. 

6.4.3 B — > VA, AA Decays to V A and AA final states have been system- 
atically studied within the framework of QCD factorization (117). The calcula- 
tions indicate that some of the tree-dominated a±p and b\p modes have sizable 
rates. For example, B(B° -► af p T ) ~ 60 x 10~ 6 and B(B° -► bfp*) « ,6(5° -► 
6^~/9 + ) ~ 32 x 10~ 6 . Likewise, the AA modes such as afa>i, afci®, a^b^ and afb® 
are expected to have branching ratios of (20 — 40) x 10~ 6 . Of these, there are only 
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a few upper limits from BABAR. The preliminary result for B(B° — > fr/ 1 /^) is 
that the branching fraction is < 1.7 x 1CT 6 (90% C.L.). The strong disagreement 
with the theoretical prediction is not understood. 

A comparison of theory with the current data on B + — > cj)Ki{\27§) + , B + — ► 
^i(1400) + (107) and B + -> a|i^* (131) seems to imply that penguin an- 
nihilation is small in penguin-dominated B — ► decays. The prediction of 
/i / (0-ftTi(127O) + ) ~ 0.44 in the absence of penguin annihilation (117) agrees well 
with the experimental result of / L (^i(1270)+) = 0.46±g;i§ (107). This indicates 
that it is the NLO correction rather than penguin annihilation that is responsible 
for pushing the longitudinal polarization fraction in B + — ► <f)K\(127Q) + down to 
the level of 0.5 . 

6.4.4 B — ► TP, TV Many charmless B decays with a tensor meson in the 
final state have been observed at B factories (see Fig. 9). Moreover, BABAR 
has measured fi in the decays B — > </>if|(1430) and B — ► wK|(1430). Contrary 
to the penguin-dominated VV modes such as (j)K* and pK*, 0X1(1430) has 
Jl ~ 0.85. Intriguingly, the B — ► wif|(1430) modes have /l consistent with 0.5. 
So far there are only two theoretical studies of the charmless B decay to a tensor 
meson, both done with the generalized factorization approach (134). Neither of 
these calculations has a prediction for fi- At the moment, the data for decays 
to tensor mesons is well ahead of the theory. 

7 BARYONIC B DECAYS 

7.1 Experimental Status 

The experimental results for 15 decays involving light baryons are summarized in 
Fig. 10. In most cases, the agreement between experiment and theory is good. 

7.2 Threshold Enhancement 

A peak near the threshold area of the dibaryon invariant mass spectrum has been 
observed in many 3-body baryonic B decays. The so-called threshold effect indi- 
cates that the B meson is preferred to decay into a baryon-antibaryon pair with 
low invariant mass accompanied by a fast recoil meson. Threshold enhancement 
was first conjectured by Hou and Soni (140). They argued that in order to have 
substantial branching fractions for baryonic B decays, one has to reduce the en- 
ergy release and at the same time allow for baryonic ingredients to be present in 
the final state. This is indeed the near-threshold effect mentioned above. 

While various theoretical ideas (140, 153) have been put forward to explain 
the low mass threshold enhancement, this effect can be understood in terms of 
a simple short-distance picture (154). In the two-body decays, one energetic qq 
pair must be emitted back to back by a hard gluon in order to produce a baryon 
and an antibaryon. Since this hard gluon is highly off mass shell, the two-body 
decay amplitude is suppressed by order of a s /q 2 . In the three-body baryonic B 
decays, a possible configuration is that the baryons in the pair B\&2 are collinear 
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Figure 10: Experimental results for decays with baryons from BABAR (blue) 
(141-144) and Belle (red) (145-150) and theoretical predictions (green) (151,152). 



and in the opposite direction from the meson. At the quark level, the quark and 
antiquark emitted from a gluon are moving in nearly the same direction. Since 
this gluon is close to its mass shell, the corresponding configuration is not subject 
to the short-distance suppression. This implies that the dibaryon pair tends to 
have a small invariant mass. 



7.3 Two-body and Three-body Decays 

None of the two-body charmless baryonic B decays have been observed so far and 
the present limit on their branching ratios has been pushed to the level of 10 -7 
for B — > pp and Ap (141, 146). The fact that three-body final states have rates 
larger than their two-body counterparts, i.e., T(B — ► BiB%M) > T(B — > B1&2) is 
due to the threshold effect discussed above. 

The study of 3-body decays is more complicated. The factorizable contribu- 
tions fall into two categories: (i) the transition process with a meson emission, 
{M\(qsq2)\0) (BxB2\(qib)\B) , and (ii) the current-induced process governed by the 
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factorizable amplitude (BiB2\(qiq2)\0) {M\(q 3 b)\B) . The interested reader is re- 
ferred to (153) for further details. 

7.4 Radiative Decay 

Naively it appears that the bremsstrahlung process will lead to T(B — ► ,61,627) ~ 
0(a em )T(B BjBa), with 

Oiem being an electromagnetic fine-structure constant, 
and hence the radiative baryonic B decay is further suppressed than the two-body 
counterpart, making its observation very difficult at the present level of sensitivity 
for B factories. However, there is an important short-distance electromagnetic 
penguin transition b — ► 57. Because of the large top quark mass, the amplitude 
of b — ► 57 is neither quark-mixing nor loop suppressed. Moreover, it is largely 
enhanced by QCD corrections. As a consequence, the short-distance contribution 
due to the electromagnetic penguin diagram dominates over the bremsstrahlung. 
The relatively large predictions of order 1 x 10 -6 (155-157) have been confirmed 
with a measurement from Belle for the decay B + — > pAj (146). 

7.5 Angular Distribution 

Measurement of angular distributions in the dibaryon rest frame will provide fur- 
ther insight to the underlying dynamics. The SD picture and the pole model 
both predict a stronger correlation of the outgoing meson with the baryon than 
the antibaryon in the decay B — ► B1B2M. This feature has been confirmed for 
B + -> ppir + and B + -> pA-f, but not for B + -> ppK+ . Both BABAR (142) and 
Belle (158) found that the K + in the latter decay prefers to be collinear with the 
p in the pp rest frame, contrary to the above expectation. This angular correla- 
tion puzzle indicates that either some long-distance effects enter and reverse the 
angular dependence or the dibaryon pair pp is produced from some intermediate 
state e.g. baryonium. 

Recently Belle has made a new measurement of the angular distribution of 
B + — > pAir + (147). Naively, it is expected that the pion has no preference for 
its correlation with the A or the proton as the dibaryon picks up energetic s and 
u quarks, respectively, from the b decay. However, the new Belle measurement 
indicates a correlation between the pion and the A. In short, the correlation 
enigma has been found in the penguin-dominated modes B + — ► ppK + and B + — ► 
pAir + and it cannot be explained by the SD b — > sg picture. This poses a great 
challenge to theorists. 

8 TIME-DEPENDENT CP VIOLATION 

The B Factories at KEK in Japan and PEP-II in California have asymmetric 
energies: the electron beam has an energy of 8 — 9 GeV while the positron beam 
is about 3 GeV. This asymmetry means that the center of mass of the T(4S I ) 
is moving so that the produced B° and B mesons do not decay at the same 
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point. Belle and BABAR take advantage of this to measure, with a precision of 
about 100 jum, the distance between the "signal" B decay and the "tagged" B 
decay. The signal final state can come from either B° or B°; the tagged B is 
either a B° or B° with the flavor determined primarily by the charge of leptons or 
kaons in the event. This measurement, together with the known properties of the 
motion of the T(45) system, allow measurement of the (signed) time difference 

5t between the B° and -B° decays. Then the time-dependent asymmetry of the 
decays to a final state / is measured: 

T(B(6t)^f)-T(B(5t)^f) 
r(B(6t)^f) + F(B(5t)^f) 



S f sm(Am5t) - Cf cos(Amdt) , (20) 



where Am is the mass difference of the two neutral B eigenstates, Sf monitors 
mixing-induced CP asymmetry and Cf measures direct CP violation (Belle uses 

Af = ~ Cf) - 

In 2001, BABAR and Belle used this technique to observe CP violation in the 
B meson system for the first time, measuring the value of S for the decay B — ► 
J/tpK° and similar b — ► ccs decays. Since there is only one quark- level process 
involved in these decays, this is known to measure sin 2/3 with an uncertainty of 
~0.001, where f3 (also called <f>\) is one of the angles of the CKM triangle (See 
Fig. 11). The ambiguity in extracting (3 is resolved with other measurements so 
that the current world average is [3 = (21 ± 1)°. In this section we discuss how 
charmless hadronic B decays can be used to measure the CKM angle a and to 
search for physics beyond the Standard Model. 




v cd v* b 

Figure 11: The CKM triangle, showing the angles a (^2), P (</>i)> and 7 (1^3). 



8.1 Measurements of at 

Measurements of S in the decay B° — ► 7r + 7r~ can be used to extract the angle 
a. However there are complications since the penguin amplitudes are not small 
("penguin pollution"), so a second weak phase is introduced. Gronau and London 
showed that an isospin analysis (159) involving also the decays B° — ► 7r°7r° and 
B + — > 7r + 7r° can be used to determine the amount of penguin pollution and hence 
extract the angle a. 

While the measurements of S and C for the decay B° — > it + it~ are fairly 
precise (5S = 0.07) (21, 160), the sensitivity to a is rather poor since the penguin 
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pollution is large in this decay. The problem can be seen from the relatively large 
branching fractions for the penguin-dominated B — > Kir decays compared with 
B — > 7T7T (Fig. 2). The situation is reversed for the VV decays where B — > ivT*p 
tends to be smaller than — > pp (Fig. 8). Thus the most sensitive measurements 
of a have come from measurements of S for the decay B° — > p + p~ (103, 161) 
and the isospin analysis involving also -B — > /j°/9° and £> + — > Some 
complications arise for these decays due to polarization and potential corrections 
to the isospin analysis due to p — uj mixing, electroweak-penguin amplitudes, and 
other isospin-breaking amplitudes. The magnitude of all of these effects is small 
compared with the present precision. While the exact numbers differ due to the 
statistical analysis used in the treatment of the results from Belle and BABAR, 
current measurements yield a value for a of about 90° with an uncertainty of 
~7°. The decay B° — > p ± ir^ L can also be used to constrain a but this is much 
less precise with present data. 

8.2 Measurements of 7 

Several methods have been employed for measurements of the CKM angle 7 that 
are theoretically clean. All of these methods use decays involving D mesons such 
as B — > DK and as such are beyond the scope of this review. The best of these 
measurements from Belle and BABAR find 7 = 76° with uncertainties of ~15°. 

Charmless hadronic B decays are also sensitive to 7 and there is a long history 
of suggested methods, all of which are imprecise with current data. The most 
successful methods have used global fits to many of the branching fraction and 
Aqy> measurements discussed in this review. A recent analysis of PV decays (17) 
finds 7 = (72 ±5)°, though this method is generally regarded to be less robust due 
to flavor SU(3) breaking and other theoretical uncertainties. Nevertheless, it is 
clear that measurements of the three angles sum to 180° within the experimental 
uncertainties of ~ 10°. 

8.3 Measurements of Penguin B° Decays 

Possible physics beyond the Standard Model has been intensively explored through 
measurements of time-dependent CP asymmetries in the penguin b — ► sqq decays 
such as B° — > (0, u, 7r°, rf)K°. In the SM, S for these decays should be nearly 
the same as the value measured for the b — ► ccs decays such as B° — > J/^K°; 
there is a small deviation at most O(0. 1) (162). In order to detect New Physics 
unambiguously in the penguin b — > sqq modes, it is of great importance to un- 
derstand SM predictions for the difference AS*/ = —f]fSf — Sb^ccs with 7]f = 1 
(— 1) for final CP-even (odd) states. The quantity Sf has been estimated in var- 
ious QCD-based approaches; the results, together with the measured values, are 
summarized in Table 1. Since Sb^ccs = 0.672 ± 0.024 (166), it is clear that AS*/ 
are predicted to be small and positive in most cases, while the experimental cen- 
tral values of AS 1 / are negative except for the K + K~Ks and KgKgKs modes. 
However the average of these measurements, 0.64 ± 0.04 (166), is less than one 
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standard deviation below the value from b — > ccs decays. 



Table 1: Mixing-induced CP violation Sf predicted in various approaches. 
The QCDF results are taken from (49, 125, 163, 164). There are two solutions 
with some of SCET predictions. The K + K~ K$ predictions and measurements 
exclude the cf> mass region. 



QCDF pQCD (28,165) SCET (26,66) 



+0.00 
-0.04 

+0.03 
-0.06 



Expt (167, 168) 



<j)K s 
loK s 

p°K s 

r,'K s 

r\K s 

7T K S 

M980)K S 
K+K-Ks 

KgKgKg 

K~ ^0„0 
JXqTT 7T 



0.75 
0.85 



0.64t^ 3 7 
0.741™° 

o.79±g : g 

7Q+ - 02 

u -' y -0.04 

n 70-1+0.001 
u - ' O1 -o.ooi 

n 7OQ+0.009 
u - ' zo -0.020 



+0.009 
-0.020 
72q +0.009 
' ^-0.020 



0.719 




0.71 ±0.01 
0.841H? 



50 +01 ° 



n 74+0.02 

u -' 4 -0.03 



0.69 
^n +0 - 05 
0.80 ± 0.02 

o.85ig:8i 

o.56i8;8i 

0.706 ± 0.008 
0.715 ±0.010 
0.69 ±0.16 
0.79 ±0.15 
0.80 ± 0.03 



0.44: 



0.18 



0.45 ± 0.24 



0.63l8il 



0.60 ± 0.07 



0.57 ±0.17 
62+ - 11 

u - oz -0.13 

0.82 ±0.07 
0.74 ±0.17 
-0.52 ±0.41 



9 CONCLUSIONS 

In this review, we have summarized branching fraction results for more than 100 
charmless i?-meson decays. Many of these decays have significant experimental 
signals. We also have shown results for CP-violating asymmetry measurements 
for nearly 50 of these decays. This represents a truly impressive body of work, 
most of which have come from the Belle and BABAR experiments in the last 
decade. 

The global features of the branching fractions and CP asymmetries of these 
charmless B decays are generally well described by the QCD-motivated theories 
such as QCDF, pQCD and SCET. The agreement between theory and experiment 
is generally satisfactory. However, there remains some unsolved puzzles: (i) Kir 
CP puzzle: it is naively expected that B + — ► K + tt° and B° — > K + ir~ have 
similar direct CP asymmetries, while they differ by 5.3 a experimentally, AAk-k = 
Acp(K + TT ) - Acp(K + tt-) = 0.148 ± 0.028; (ii) the abnormally large B -► rfK 
rates: while the qualitative picture of the enormously large rate of B — > rj'K 
over B — > r\K is understood, a precise quantitative prediction is still lacking; 
(iii) branching fractions of B° — > 7r°7r°, p°7r°: theory usually predicts too small 
rates for them and a reasonable one for B° — > p°p°; (iv) the polarization puzzle 
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in penguin-dominated B — > decays: the transverse polarization fraction is 
not as small as naively anticipated; (v) mixing-induced CP asymmetries: the 
predicted values of the effective sin 2(5 for most of the b — ► sss induced decays 
are above the one obtained from B — > J/tpK , whereas experimentally the value 
of sin 2/3 from the bulk of the decay modes is systematically below that of B —>■ 
J/ipK°; (vi) the angular correlation enigma in three-body baryonic B decays: the 
short-distance picture of b — > sg cannot explain or accommodate the observed 
angular distributions in penguin-dominated decays such as B + — ► pK + , pAn + 
and B° -► AAK°. 

The aforementioned puzzles pose a great challenge to the B-physics community 
and their solutions need efforts from both theorists and experimentalists. Either 
these enigmas can be resolved within the framework of the standard model pro- 
vided that the hadronic matrix elements are under fair control or new physics 
effects already manifest themselves in some of the puzzles. 
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